Strongest chemical weathering in response to the coldest period in Guyuan, Ningxia, China, during 14-11 Ma

Moisture evolution in Central Asia including Northwest China shows less similarity with its surroundings and attracts a growing number of studies. In this study, a well-dated thick lacustrine sequence is chosen in Northwest China and detailed geochemical analysis is conducted during the Middle Miocene Climate Transition (MMCT, 14–11 Ma). The multi-proxy records (Na2O/Al2O3, CIA, Rb/Sr) revealed that chemical weathering was the strongest during 11.85–11 Ma, the coldest period in 14–11 Ma as evidenced by the global deep-sea oxygen isotope records. Accordingly, we conclude that global climate cooled during MMCT and reached the coldest during 11.85–11 Ma. Thus, the westerly circulation became the strongest during this period, which brought more water vapor to Northwest China and the chemical weathering was significantly improved. On the other hand, the significant decrease in temperature led to the marked weakening of evapotranspiration, and thus the effective humidity was relatively increased. Both aspects contribute greatly to the significant enhancement of chemical weathering in eastern Central Asia. This weathering history of the sediments in the northeastern Tibetan Plateau is of great scientific significance to understanding tectonism and climate change in Asia during MMCT.


Introduction
During the Mid-Miocene Climate Transition (MMCT) period, the earth's climate experienced long-term cooling, the sea level dropped significantly and the polar ice sheet increased sharply [1][2][3][4]. Continental records in Asia show a good response to global cooling during 14-11 Ma, and indicate that global cooling led to not only weakening of the summer monsoon and declining of vegetation cover [5][6][7][8], but also strengthening of the winter monsoon as evidenced by coarsening of eolian dust particles [9][10][11].
Noticeably, recent studies indicate that moisture evolution in Central Asia including Northwest China shows less similarity with its surroundings, mainly expressed as moisture increase with global cooling [11,12]. The topographic change resulting from the Tibetan Plateau uplift was proposed to have strongly influenced the moisture patterns in Central Asia during MMCT  [11]. The decrease in sedimentary leaf wax δD n-alk between 15 Ma to 10. 4 Ma was used to estimate the gain in elevation ranges between 1.6 km and 2.5 km with a mean of 2.1 km [12]. However, analysis of δ 18 O data from 2750 sedimentary carbonate samples across Asia suggests that a long-standing topographic feature has continuously blocked southerly moisture and subsequent progressive uplift of the Tibetan Plateau has had little impact on the Central Asian climate [13]. Therefore, it remains disputed whether the moisture increase with global cooling during MMCT is attributed to tectonic uplift or the westerlies strengthening. It is well known that chemical weathering intensity is largely controlled by temperature and precipitation [14,15]. High precipitation and/or warm temperatures can enhance the chemical weathering intensity, whereas either low temperature or decreased precipitation can decrease the chemical weathering intensity [14]. In this study, an analysis of chemical weathering was conducted on the fluviolacustrine sediments to examine the role of westerlies strengthening on the moisture increase in Central Asia during MMCT.
In this study, a well-dated thick lacustrine sequence is chosen in the eastern Liupan Mountains [16]. Given that the Mid-Miocene Climatic Optimum (MMCO) occurred at 16-14 Ma in the Longzhong Basin [17], the detailed geochemical analysis is conducted from 14 to 11 Ma. Our aim is to identify changes in sediment origin and weathering history for the sediments in the northeastern Tibetan Plateau during sedimentation. It is of great scientific significance to understand tectonism and climate change in Asia during the late Middle Miocene.

Geological and geographical setting
The Sikouzi area is situated on the east side of the Liupan Mountains, approximately 40 km northwest of the town of Guyuan, Ningxia Province, and has a mean elevation of about 1550 m a.s.l (Fig 1). The north of this area is surrounded by the Tengger, Wulanbu, and Mu Us Deserts. Influenced by the East-Asian summer monsoon, climatic conditions of the Sikouzi area at present are temperate and characterized by relatively hot, humid summers and cold winters. For the past 30 years, the mean annual temperature (MAT) is 6.2˚C, with a July average of 18.9˚C and a January average of -8.3˚C. Mean annual precipitation (MAP) is 478 mm, and over 60% of the precipitation falls in July-September with a peak mean rainfall of 109.1mm in August [5,10,18]. Mean annual latent evaporation reaches 1772 mm, which is 3.7 times the annual precipitation. The East Asian winter monsoon, primarily controlled by the Siberian High, drives strong northwestern winds below 2000 m altitude mainly from December to April in this area [5,10,18].
The Sikouzi area is a transitional region between the northeastern margin of the Tibetan Plateau and the Ordos block, which has been relatively stable since the Cenozoic era [19,20]. In the southwest is the Haiyuan-Liupanshan arc fault zone with strong tectonic deformation in the late Cenozoic [20][21][22], and the Xiangshan-Tianjingshan fault is the boundary in the northeast.
In the study area, reddish and brownish fluviolacustrine sediments, with a thickness of 2880 m, are well exposed along a 5265 m stretch of the Qingshui River. The cross-section map and several photographs of the typical outcrops are presented in Jiang et al. [16]. The lacustrine succession of the Sikouzi area, resting pseudoconformably on Eocene sandstones, strikes NNW and dips ENE with an inclination of 18˚-64˚. The whole succession includes an anticline and a syncline and is covered unconformably by last glacial loess deposits [16].
Previous studies have suggested that the deformation of the Late Cenozoic strata in this region was caused by the Haiyuan left-slip fault, the recent activity of which began at about 0.2 Ma ago [20,21,23]. The folding caused some difficulty with stratigraphic correlation in the lower part of the section. Fortunately, thick marker beds of white sands in the lower section are distinguishable and the stratigraphic correlation was accurately completed using these marker beds. The lithology of 14-11 Ma is composed of lacustrine shore and deltaic sediments [24], which mainly consist of alterations of greyish and reddish sandstone layers. Thin gypsum layers occur occasionally within the sediments.

Methods
The exact latitude and longitude (GPS data) of the beginning and end point of the Sikouzi fluviolacustrine section is 36˚18 0 N, 106˚02 0 E and 36˚16 0 N, 105˚59 0 E, respectively. In the field, samples were carried out at intervals of 10-30 m along the Sikouzi section, and there is no need for any permits for field or sample access. Sampling sites were selected where possible in mudstone, silty mudstone, muddy siltstone and siltstone. At each site the surface weathered material was removed and a fresh sample was taken. The samples were stored in the State Key Laboratory of Earthquake Dynamics, Institute of Geology, China Earthquake Administration. Age control for all the samples is derived from a detailed palaeomagnetic record and biostratigraphic data [16]. In this study, thirty samples were selected strictly at an interval of 0.1 Ma spanning 14-11 Ma after age was obtained.
Major and trace element concentrations of bulk samples were determined at the Institute of Geophysical and Geochemical Exploration, Chinese Academy of Geological Sciences following the method of previous studies [25,26]. The samples were fully powdered using an agate mortar, and pretreated to remove organic substances and carbonate components. 10% Hydrogen peroxide was added to remove organic substances. 1 M acetic acid (HOAc) was added to remove carbonate components. Each step was stirred sufficiently, respectively. The solid residues were cleaned by centrifugation three times in pure water. The samples were dried in an oven at 40˚C and fully powdered. Approximately 4 g of ground sample was weighed, transferred with boric acid into the center of a column apparatus, and pressurized to 30 t/m2 for 20 s using the ultra-high pressure sample (UHPS) preparation system. The compressed samples, approximately 4 cm in diameter and 8 mm in thickness, were analyzed by PW4400/40 X-ray fluorescence spectrometer. About 0.7 g of ground sample was weighed in a crucible, and put in a muffle furnace to calculate the loss on ignition (LOI).
The most significant feature of this record is that the abundance of CaO reaches 3.51% at about 11.85 Ma, which is much higher than the average percentage of 0.58% of the whole sequence during 14-11 Ma. Meanwhile, the percentage of Al 2 O 3 (6.34%), K 2 O (2.09%) and MnO (0.02%) are the lowest values in the whole record. At the same time, it is noticeable that the Na 2 O/Al 2 O 3 ratio reaches the maximum value (0.37), CIA arrives at the minimum value (36.6), Rb/Sr runs to the lowest value (0.01) of the whole record. These indicate extremely weak chemical weathering intensity in the study area.

Eolian origin of the Sikouzi fluviolacustrine fine sediments
It can be seen that all variations for these major elements fluctuate within narrow ranges (Fig 2,  Table 1). Such a pattern of the major elements is similar to those observed in the loess-paleosol deposits [28] in the CLP (Fig 3). This good exponential linear correlation suggests that the fine-  (Fig 4). This means that the finegrained sediments in the study area are not only windblown origin, but also possibly have similar dust sources to those of the loess-paleosol sequence in the CLP.

Age model assessment of the late Miocene in the Sikouzi section
Our detailed magnetostratigraphy research was published in 2007 [16]. Subsequently, Wang [22] and Lin [33-35] verified the age model again. After careful comparison, it can be found that the differences in age models for the entire section are very small, and these are even smaller during the late Middle Miocene.
There are several constraints on the comparison of the late Middle Miocene age models. Chron C5n.2n spans about 1 Ma and appears at similar depths in different age models [16,22,[33][34][35]. Tectonic activities occurred in the Sikouzi, Dahonggou and Linxia basin at the almost same time of~12 Ma [5,36,37] though with a more significant intensity at~8 Ma around the Liupan Mountains [38 -43]. The end-member analysis of all 3499 grain-size samples in Sikouzi lacustrine sediments indicates that the varying trend of three end-members can be successively correlated in seven stages with the integrated benthic δ 18 O record [44]. The magnetostratigraphic record of the Sikouzi section was not only identified the all polarity chron and polar sub-chron but also constrained by the hipparion fossils [16]. Furthermore, the enhanced chemical weathering intensity indicated by the multi-proxy records (Na 2 O/Al 2 O 3 , CIA, Rb/Sr) corresponded well in timing with the cooling period revealed by the global deepsea oxygen isotope records [45] during 14-11 Ma, and there are 9 events of intensified chemical weathering correspond well to the cold periods evidenced by the deep-sea oxygen isotope records [45] within the dating error in this period. (as described below). Based on the above considerations, we still use the previous age model [16].

Geochemical proxies and their environmental implications
The chemical index of alteration (CIA) is often used as a good indicator of the degree of weathering of sediments. It is calculated in molecular proportions as follows: CIA = [Al 2 O 3 /(Al 2 O 3 + CaO � + Na 2 O + K 2 O] × 100, where the CaO � is the amount of CaO in the silicate minerals [27]. In the process of chemical weathering, Ca, Na and K in feldspar (main silicate) are easy to leach compared with Al. With the enhancement of weathering, the contents of Ca, Na and K in weathering products is reduced, while the contents of Al remained relatively unchanged. Therefore, an increase in the CIA value indicates the enhancement in chemical weathering, and the amelioration of climate [46]. Similarly, the ratio of Na 2 O/Al 2 O 3 can be used as an indicator of chemical weathering [46,47], as the mother rock undergoes chemical weathering, in which the unstable Na + is leached and the stable cation Al 3+ is basically unchanged and relatively enriched. Therefore, the smaller Na 2 O/Al 2 O 3 ratio, the stronger the chemical weathering, and vice versa.
The Rb/Sr ratio is an important indicator of weathering intensity in both lacustrine and loess sediments [48]. In loess sediments, the loss of Sr content is mainly controlled by carbonate loss, Rb may be lost during the transformation of clastic mica and other potassium-bearing minerals into illite and some clay minerals, but the clay minerals can also contain Rb in their structures, so the loss of Rb in paleosoils is much smaller than that of Sr [49]. Therefore, the ratio of Rb/Sr can reflect the strength of chemical weathering of loess strata, high Rb/Sr ratio indicates strong chemical weathering and warm and humid climate while low Rb/Sr ratio indicates weak chemical weathering and cold and dry climate [46,49,50]. Similarly, the high Rb/Sr ratio of lake sediments can reasonably explain the strong chemical weathering in the lake catchment [49,51].

Possible tectonic activity at around 11.85 Ma.
The Sikouzi chemical weathering record can be divided into two parts by 11.85 Ma (Fig 5). There is a significant difference between the upper part  and the lower part . From the lower to the upper part, the mean Na 2 O/Al 2 O 3 ratio decreases from 0.18 to 0.11, CIA increases from 66.7 to 69.9 on average, the mean Rb/Sr ratio increases from 0.8 to 1.59. Comparably, from the Sikouzi sequence, the mean L � value increase from 53.5 to 57.6, the mean a � value decrease from 13.1 to 10.4. What's more, deep-sea oxygen isotopes increase from 2.32 to 2.57 on average [45].
Intriguingly, the significant change at 11.85 Ma in the study area however has no evident response from the deep-sea oxygen isotope record [45]. On the other hand, several previous studies reveal that tectonic activity occurred widely in the northeastern Tibetan Plateau at the almost same time of around 12 Ma. Carbonate oxygen isotope analysis of lake and river sediments in the Linxia basin shows that, till 12 Ma, the Tibetan plateau has risen high enough to  12 Ma [37,52]. More records revealing tectonic activity at around 12 Ma in and around the Tibetan Plateau are summarized by Ma and Jiang [53] and in their reference. Accordingly, we consider that a tectonic event occurred around 11.85 Ma in the study area, which can be well correlated and linked with the tectonic change at around 12  Correspondingly, the L � value (53.29-63.27, average 57.59) in 11.85-11 Ma reached the maximum of the whole sequence, indicating the maximum carbonate content [18]. The a � value (8.09-12.06, average 10.42) in 11.85-11 Ma reached the lowest of the Sikouzi sequence, revealing the lowest temperature [18]. These results can be compared with the maximum value of the deep-sea oxygen isotope record [45] during 11.85-11 Ma in the Middle Miocene Climatic Transition (MMCT, . This indicates that the chemical weathering in the study area is the strongest when global ice volume reaches the maximum and the temperature reaches the minimum during MMCT.  Fig 5). Obviously, the chemical weathering intensity weakens relative to that in 11.85-11 Ma. Noticeably, there are 9 events of intensified chemical weathering in this period, and their intensities are close to that in 11.85-11 Ma. Intriguingly, these 9 events correspond well to the cold periods evidenced by the deepsea oxygen isotope records within the dating error [45].

Possible mechanism of increased chemical weathering in cold period in eastern Central Asia
High precipitation and/or warm temperatures can enhance the chemical weathering intensity [14]. Considering that 11.85-11 Ma is the coldest period in the whole sequence but with the strongest chemical weathering, we believe that moisture is the most abundant in this period. Generally, there are three hypotheses for the abnormal increase of moisture in Central Asia during global cooling, retreat of the Paratethys, global cooling, and uplift of the Tibetan Plateau.
The proto-Paratethys Sea covered a vast area in Central Asia during the late Eocene and significantly influenced regional climate by providing an immediate source of water vapor [54][55][56]. However, recent studies show that there were three obvious transgression/regression cycles of the proto-Paratethys Sea [57,58]. The first (from~59-57 Ma to~53-52 Ma) and the second (from~47-46 Ma to~41-40 Ma) incursion have been poorly constrained, while the third incursion has been precisely dated to 39.7-36. 7 Ma. This implies that the final retreat of the Paratethys Sea occurred at 36. 7 Ma, and since then, the climate of eastern Central Asia was little affected by the retreat of Paratethys Sea [59].
Uplift of the Tibetan Plateau in the late Cenozoic has generally been believed to have played a significant role in strengthening the Asian monsoon through modulating the atmospheric circulation and its barrier effect to southern-sourced moisture [18]. However, the timing of the uplift of the Tibetan Plateau is controversial. Thermochronometry studies indicate that rapid exhumation of the northeastern Tibetan Plateau started at~10 Ma in the Qilian Shan, at 8 Ma in the liupan Shan [38][39][40][41][42], at 10-8 Ma along the Kunlun and Haiyuan faults [60] and at 13 Ma in the Jishi Shan [61]. This uplift timing slices are apparently later than 14-11 Ma. In addition, a compilation of δ 18 O records across Asia suggest that progressive uplift of the Tibetan Plateau has had little impact on Central Asian climate [13].
Global cooling exerts a major effect on changes in water vapor in Central Asia in two aspects. On one hand, global cooling leads to weakening of evaporation and transpiration and consequently increases water vapor content relatively in the study area. On the other hand, global cooling leads to enhancement of westerlies circulation, which transports more Atlantic water vapor to Central Asia. Therefore, we believe that during 14-11 Ma, global cooling led to enhancement of the westerlies circulation, which brought more water vapor into the study area and resulted in the increase of chemical weathering intensity as during 50.5-37. 6 Ma in eastern China [47]. This recognation is supported by the improvement of vegetation and climate conditions in the Xunhua Basin during 12.5-8.0 Ma, which was attributed to the decrease in evaporation rates caused by continuous global cooling [62].

Conclusion
The multi-proxy (Na 2 O/Al 2 O 3 , CIA, Rb/Sr) of the fluviolacustrine sediments at Guyuan, Ningxia, China during MMCT  revealed that chemical weathering was the strongest during 11.85-11 Ma. The global deep-sea oxygen isotope records show that temperature during 11.85-11 Ma was the lowest in 14-11 Ma. Accordingly, we conclude that global climate cooled during MMCT and reached the coldest during 11.85-11 Ma. Thus, the Northern Hemisphere climate gradient became the largest and the westerly circulation became the strongest during this period, which brought more water vapor to Northwest China and the chemical weathering was significantly improved. On the other hand, significant decrease in temperature led to the marked weakening of evapotranspiration, and thus the effective humidity was significantly increased. Both aspects contribute greatly to the significant enhancement of chemical weathering in eastern Central Asia.